Formal Solution of the Fourth Order Killing equations for Stationary Axisymmetric 

Vacuum Spacetimes. 
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An analytic understanding of the geodesic structure around non-Kerr spacetimes will result in a 
powerful tool that could make the mapping of spacetime around massive quiescent compact objects 
possible. To this end, I present an analytic closed form expression for the components of a the 
fourth order Killing tensor for Stationary Axisymmetric Vacuum (SAV) Spacetimes. It is as yet 
unclear what subset of SAV spacetimes admit this solution. The solution is written in terms of 
an integral expression involving the metric functions and two specific Greens functions. A second 
integral expression has to vanish in order for the solution to be exact. In the event that the second 
integral does not vanish it is likely that the best fourth order approximation to the invariant has 
been found. This solution can be viewed as a generalized Carter constant providing an explicit 
expression for the fourth invariant, in addition to the energy, azimuthal angular momentum and 
rest mass, associated with geodesic motion in SAV spacetimes, be it exact or approximate. I further 
comment on the application of this result for the founding of a general algorithm for mapping the 
spacetime around compact objects using gravitational wave observatories. 



I. INTRODUCTION 

During an extreme mass ratio in-spiral (EMRI), a 
small probe object such as a neutron star spirals into 
a much more massive compact object commonly held to 
be a black hole while broadcasting a distinctive gravi- 
tational wave (GW) signal. The detection of these GW 
signals using detectors such as LIGO and LISA could in 
principle allow us to read off the structure of the massive 
compact objects 0, HI, much as the head of a gramo- 
phone player reads the tune off a long playing record. 
The feat of so mapping spacetime has been shown to be 
possible in principle 0] , however an explicit algorithm for 
doing so in practice for generic objects is still lacking 
The GW signal emitted during an inspiral depends on 
the trajectory of the probe, as well as the propagation of 
the GW field away from the compact object. This trajec- 
tory can be approximated as a flow through a sequence of 
geodesies and one of the main stumbling blocks hindering 
the development of an algorithm for mapping spacetime 
is finding an explicit description of the geodesic structure 
of the spacetime surrounding the object [4[. 

Our knowledge of the geodesic structure of spacetimes 
to date encompasses only a subclass of Petrov type D 
spacetimes found to be separable by Carter 0, Q . This 
class of spacetimes, which includes the Kerr metric de- 
scribing the gravitational field of a rotating black hole, 
has many special properties @, 0, H, 0] which currently 
undergird our approach to generating EMRI waveforms 
[Tol | . A key feature of the Carter spacetimes is that 
they admit a second order Killing tensor T( 
obeys the Killing equations (KE), Ti 



(Q1Q2) 



which 



(a 1 a 2 ;a 3 ) 



0. As a 



result, for a particle with momentum p a , the quantity 
Q2 = T( aiC " 2 ' > p ai p a2 remains constant along its trajec- 
tory. This Carter constant Q2, in addition to the con- 
stants related to conservation of rest mass \i, energy E 
and azimuthal angular momentum L z uniquely deter- 
mine the geodesic structure of the Carter spacetimes. 



General SAV spacetimes are of Petrov type I and do 
not admit second order Killing tensors [ijj. However, 
numerical integrations [HI, EH performed on SAV space- 
times, exhibited an orbital crossing structure indicative 
of the existence of a fourth order orbital invariant Q4 
quartic in momenta [l3T |. Such an constant of motion 
can be expressed as 

W4 — 1 Pa\Pa 2 Pa 3 Paa.i ) 

where the Killing tensor T obeys the fourth order KE 



± (aia2a 3 a4;a^) u - 
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It should be noted that third order Killing tensors need 
not be considered as discussed in [IH . When written out 
in full the fourth order KE (|2|) represent 56 coupled dif- 
ferential equations, governing the 35 independent com- 
ponents associated with a totally symmetric tensor T. 
Possibly for this reason, no previous analytic work relat- 
ing to the solution of fourth order Killing tensors appears 
to exist. 

This paper provides a non-trivial closed form solution 
to the fourth order KE ||2J| in SAV spacetimes. It is the 
culmination of a serious of papers [4l.l6l.llll.ll3l| investigat- 
ing the requirements for mapping spacetime and finding 
a suitable geodesic description for SAV spacetimes that 
could make such an endeavor tractable. The formalism 
and variables used in this derivation were originally de- 
signed to produce an algebraic check for the existence of 
fourth order Killing tensors and are discussed more fully 
111 It came as a considerable surprise that one could 
actually write down the closed form solution presented 
here. 

The analysis is valid for all SAV spacetimes with two 
commuting Killing vectors, dt and d$, represented by 
means of the Lewis-Papapetrou metric 
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The metric functions are entirely determined by a com- 
plex Ernst potential £ [ljj , the real part of which is 



terms of the metric functions as follows: 



2V 



Line integrals of £ determine the functions 



7 and u. The function R is any harmonic function obey- 
ing the equation R zz + R pp = and represents a residual 
coordinate freedom in the metric. 

In Sec. |n] a "symmetric" formulation of the fourth 
order KE derived in [ll[ is given. Some of the features 
which make explicit solution of these equations possible 
are highlighted. A constructive derivation of the sec- 
ond order Killing tensor components associated with the 
Carter spacetimes, given in [6j, indicated that it is possi- 
ble to write down a formal solution to the equations be- 
fore determining the explicit SAV metric functions that 
admit such a solution. A similar feat for the fourth order 
KE is performed in Sec. IIIIl The resulting solution is 
given in the form of an integral expression involving the 
metric functions and appropriately chosen Green's func- 
tions. In addition the requirements for the existence of 
the solution are also given in integral form. In conclu- 
sion the implications of the given solution for a general 
algorithm for mapping spacetime are discussed in Sec lIVI 



II. FOURTH ORDER KILLING EQUATIONS 

The KE for SAV spacetimes were carefully analyzed in 
fj~fl ] and an alternative "symmetric" formulation of the 
fourth order KE derived. This formulation can be viewed 
as seeking a solution to four "interlocking" fourth order 
Killing tensor problems for a two-manifold and is entirely 
equivalent to Eq. in the SAV case. The explicit linear 
transformation between the non-zero fourth order Killing 
tensor components T < ' aia2a3a4 " > and the variables, -P<i : j> 
used in this section are given in Sec. VIII of [1JJ. (The 
transformation constitutes roughly a page.) The formu- 
lation of the KE in terms of P<i:j> and the resulting 
structure and simplicity is pivotal in making the prob- 
lem of solving the large set of KE analytically tractable. 

The variable names P<i : j> are based on the derivative 
structure imposed by the KE and can be divided into two 
sets. In the first set, the indices j ' — 3 and j ' = 4 indicate 
that the KE fix derivatives with respect to £ = l/2(p+iz) 
and C = l/2(p — iz) respectively. The index i labels 
different variables with a fixed derivative structure. The 
second set with index j = indicates that the KE fully 
determine the gradient of the variable. There are nine of 
these variables, P<ko> with i G {—4 1, 1 ■ • ■ 5}. 

The KE governing the derivatives of the 10 variables 
of the first set are given by the linear system: 



P<z:3>, C = /,-/' 5:3 ~ ^ P <5:3>X * e {I"'*}, 
P < t A>X = -fi,C P <5-A> ~ \frP<5:4>, C »e{l-4}, 
P<5:3>,C = P <5A>X = °- ( 3 ) 



The functions fi entering these equations are defined in 



fi = = V, f 2 = 



h 
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h = 
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Note that P<5:4> in Eq. ([3]) is an analytic function of £ 
and indicates a gauge freedom still present in the metric. 
Without loss of generality, as was done in the example 
calculation for second order Killing tensors Q we can set 
P<5:4> = -P<5:3> = 1 further simplifying Eqs. ([3]) to 



P. 



<i:3>,C 



p 



<i:4>X 
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(5) 



To keep the remaining group of KE as succinct as pos- 
sible and to elucidate the structure found in them more 
fully, define the operators fp and fp to be 

f P (i,j) = -2(P <j:3 >/ f ), ? + /./';:, 

JpJlJ) = -2(P <j:4> / i ), C + fiP<yA>X (6) 

The equations governing the second set of real Killing 
variables P<j : o>, whose gradients are fully determined, 
are given below 



P<-1:0>, C ^gMM), 



P<-2:0>, C =/p(l,2) + /p(2,l), 

P<-3:0>, C =MM) + /p(3,i), 

P<-4:0>,C =/p(l,4) + /p(4,l), 

P<i:0>, C =4/p(2,2), 



p. 



<2:0>,C 



2(/p(2,3) + M3,2)), 



P<3:0>,c - - (/P(2, 4) + /p(4, 2)) + 2/p(3, 3), 
P<4:0>,c =4(/p(3,4) + /p(4,3)), 
P<5:0>,c=4/p(4,4), (7) 

The complex conjugates of Eqs. ([7]) constitute the re- 
mainder of the KE that have to be satisfied and are de- 
noted by Eqs. CC® 

A non-trivial solution of Eqs. (©, (7} and CC(J7J>, sub- 
stituted into the transformation from P<j : j> to T given 
in Sec. VIII of [ll( yields an explicit expression for an 
invariant of geodesic motion of the form (fi}. I n Sec. IIIIl 
such a solution is proposed. This solution includes the ex- 
isting class of Carter spacetimes or spacetimes admitting 
reducible fourth order Killing tensors as well as extending 
our knowledge to irreducible fourth order Killing tensors. 



III. FORMAL SOLUTION OF KILLING 
EQUATIONS VIA A FOURIER METHOD 

The KE as expressed in Sec. [TT| are governed by only 
two distinct differential operators. This feature is key to 
finding an analytic solution. These operators include the 



3 



Cauchy-Riemann operator = d p — id z used in Eq. ([5]) 
and on the left-hand side of Eq. (0) and the bilinear op- 
erator fp(i,j) defined in Eq. ([6]) that appears in various 
linear combinations in Eq. ([7]). 

Writing down the solution to the KE as given in SeclTTl 
is a three step process. We first solve for the complex 
valued functions P<i : 3> and P<iA> that obey Eq. (5]). 
These results are substituted into Eqs. ([7]), which are 
in turn also solved, yielding an integral expression for 
P<i-.o>- The solution satisfies Eqs. (0) and ([7J and need 
not be a real valued function. However if it is actually 
a solution of the overdetermined set of KE, Eqs. CC([7]) 
also have to be satisfied and this provides the constraint 
that the imaginary part of P<i : o> must vanish. 

There are several analytic solution methods [TBI, HH, E3] 
that can be used to solve Eq. (0, each yielding in- 
sight about existence, uniqueness and general behavior of 
the solution. The approach followed here is the Fourier 
method, which has the advantage of making the next step 
of solving ([7]) easily tractable. In this method the bound- 
ary conditions are implicit, the angular variable z is taken 
to have a range < z < 7r with z = n/2 indicating the 
equatorial plane and p corresponding to an asymptoti- 
cally radial variable extending to infinity. We thus take 
a Fourier series in z and a Fourier integral transform in 
p and define this operation to be the Fourier transform, 
FT. The FT of a function h(p,z) is denoted by H(£,T]) 
and defined by 



W(£,T7) = ~ 



oo JO 



-i2-qz -2i£p 



h(p, z) dp dz, (8) 



is obtained if we exchange summation and integration 
and write the solution in terms of a Greens function, 
G{p, z, p, z). In the case of -P<i : 3> this yields 



2 /"OO />7T 

P<i-.3>(P,Z) = - / / fi{P,z)G(p,Z,p,z) 
71 J-ooJo 



dp dz, 



(11) 



with the Greens function given by 



G(p,z,p,z) = -- e l2 ^ z -^ I e^P-^dt 
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(2<£ + 277) 



(12) 



Making use of the identities 



= 2i?(Ap) 



<% = ttS(A P ), J2 e 



i2fj(Az) 



irS(Az), 
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* J -00 2(i£ + n) 



di = H{p)e~ 2n ~ p if n > and 



fj=i 



J2vx 2n = [x- -) if < 1 



(13) 



where H(x) is the Heaviside step function and 5(x) the 
Dirac delta function, the expression for the Greens func- 
tion, given in Eq. (|12[) . can be reduced to 



where £ is a continuous variable and r\ takes on integer 
values. The original function is restored by taking the 
inverse transform, IFT, which can be expressed as 



7T 



E 



J2rjz 2i£p 



By applying the FT to Eq. 
— ^+ i2£ — 2ri and iz 

m - 2»?) 



e**P H&vW. (9) 

([5]) and recalling that 
+ 2r\ we obtain 



V. 



<i:3> 



m + 2r,) 



<i:4> 



(2zg + 2 V ) 
(2i£ - 2r,) 



(10) 



where Ti represents the FT of the functions fi defined 
in Eq. (3). Taking the IFT's of (10]) yields an explicit 
expression for P<i : j>. Greater insight into the solution 



G(p, z, p. z) = -ttS(z - z)5(p — p) + 



,sinh(2C-2C), 



(14) 



Given expressions (jTTJ) and (fT5]) for the components 
-P<i:3>(/°> z ) we can proceed to the next step of solving 
Eqs. (J7J that are governed by the operator fp{i,j) de- 
fined in Eq. ([6]). All Eqs. (UJ can be written in the form 



P<k,o>,( = ^2C„fp(i ni j n ), 



(15) 



where the integers n, i n and j n and the rational constants 
C n depend on the specific equation being evaluated. Tak- 
ing the FT of Eq. (U]) yields 



z)e^ l2riZ e- 2i ^ 



^V- ! • i>- *) , S + j- (G{p, z, p, £)) \ 1 f^z). 



{2i£ + 2rf) &^ 



m + 2*7) 



(16) 
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We obtain an expression for P<k : o> by taking the I FT of Eq. (|16|) and reordering integration and summation to give 



oo JO 



dp dz f in (p, z)f jn {p, z)K(p, z,p,z, p, z) (17) 



where the "Greens" function K was calculated to be an exact Killing tensor), the real part of -P<;.o> along 



K(p, z, p, z, p, z) = 2G{p, z, p, z)G(p, z, p, z) 

+ 4 (G(p, z, p, z)) ( \=— j 

(18) 

using Eq. (fT3| . the expression for G, Eq (fT4| . and the 
identity 



-j oo 



J2ij(z-z) 



oo 2i£(p-p) 



(2^ + 2rj) 



17^ = 



1 



1 - e 4 C- 4 C 

(19) 



The expressions for the variables P<j,j> given in 
Eqs. (fTTj) and (jTTJ) with the Greens functions defined 
by Eqs. (fTT)) and (fT5|) constitute a valid solution to the 
Eqs. ||SJ and ([7]). The full set of KE is however over- 
determined and requires that the additional Eqs. CO© 
also be satisfied. This is synonymous with requiring that 
the imaginary part of Eq. (|17p is zero. The number of 
spacetimes that have this property has yet to be deter- 
mined. 



IV. CONCLUDING REMARKS 



with P<jj>, j G {3,4} given in Eq. (fTTj) should yield 
the best fourth order approximation to the invariant and 
thus provide an analytic handle for describing geodesies 
in all SAV spacetimes . 

In H a program for mapping spacetime around com- 
pact objects with arbitrary multipole moments was set 
forth. One of the key stepping stones toward attaining 
an implementable algorithm was the ability to describe 
a so-called "geodesic equivalence class" between space- 
times that allowed us to explore features of the orbital 
frequencies and the drift of these frequencies without a- 
priori knowing what the spacetime was [4|. The Green's 
function formulation of the solution to the fourth order 
KE given here allows us to do just that. It provides a 
rigorous mathematical substructure on which the sub- 
sequent algorithmic steps of, translating the GW field 
away from the compact object, and the detection phase 
can be built without knowing the specific details of the 
spacetime under consideration. 

The solution given here makes explicit that missing 
analytic link, between the geodesic behavior of a particle 
within a SAV spacetime and the structure of the space- 
time itself. The full exploration of its implications, and 
further development of the ideas set out in Q for map- 
ping spacetime will be the subject of future work . 



This paper provides a formal solution to the fourth 
order KE for SAV spacetimes. The word formal is used 
in the sense that a lot of work regarding the evaluations 
of the integrals presented here still has to be performed 
before the scope of validity of the solution is determined. 
To date the expressions given here have been successfully 
used to predict the Poincare map for the Zipoy-Voorhees 
(ZV) metric [Hj]. The analytic results agree with well 
over 500 numerical orbits integrated so far. 

In the event that the imaginary part of P<i,o> given 
in (|17p is non-zero, (and thus the solution given is not 
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